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Abstract 

Under an appropriate change of the perturbation variable Lifshitz- 
Khalatnikov propagation equations for the scalar perturbation reduce to 
d'Alembert equation. The change of variables is based on the Darboux 
transform. 

1 Introduction 

The gauge-invariant perturbation theories efficiently eliminate non- 
physical perturbations, yet they provide the propagation equations in noncanon- 
ical form — different for each theory. In this case, pure artefacts of the choice 
of the reference system, gauge or the perturbation variable, are likely to be 
confused with new dynamical phenomena yet unknown in the laboratory-scale 
physics. This particularly concerns the large scale limit, where some "non- 
oscillatory behaviour" outside the particle horizon is commonly expected. 

In this paper we show that the definition of the gauge- invariant variables still 
possesses a freedom to chose the canonical variables, i.e. variables which satisfy 
d'Alembert equation in its standard form. To show that we have used the classi- 
cal Lifshitz-Khalatnikov formalism, which after being appropriately extended is 
equivalent to other "manifestly" gauge invariant descriptions. A systematic con- 
struction of the gauge-invariant canonical variables is provided. Consequently, 
the scalar perturbation propagates like a massless scalar field in the Roberson- 
Walker space-time, and therefore there is a close analogy between the pertur- 
bation theory (the acoustics of the expanding universes) and the field theory in 
the curved space-time. Propagation of the sound waves in the early universe 
may be also considered as an example of the acoustic geometry in the sense of 
Unruh [T]. 

The paper is organized as follows: In the second Section we remind the 
method of elimination of spurious modes in the Lifshitz-Khalatnikov theory 0. 
This Section is a generalization of the procedure given in to the case of 
arbitrary space curvature and arbitrary density-dependent pressure p = p(e). 



In Section [3] we present the general method of construction of gauge-invariant 
quantities in the synchronous system of reference. We show that the procedure 
is not unique, and therefore, there exist vast classes of physically relevant gauge- 
invariant variables. The remaining freedom to choose between them may be used 
to better describe the perturbation dynamical properties, or to construct more 
adequate observables. A specific choice depends on the researcher motivation, 
and on the character of the problem to be solved. Finally in Section0]we adopt 
the previously defined methods to obtain a canonical form of the perturbation 
equations — the d'Alembert equation. 



2 Synchronous system of reference — techniques 
of reduction of the gauge freedom 

Consider the Robertson- Walker universe of arbitrary space curvature [K — 
-1,0, +1) 



g^u = a 2 (r?) diag 



K 



K 



with the hydrodynamic energy-momentum tensor 

= (e + p)u^u" +pg fiu 



(1) 



(2) 



and the barotropic equation of state p = p(e). We investigate small perturba- 
tions Sg^ of the metric (y. We limit ourselves to scalar (density) perturbations, 
therefore, 8g^ v is defined by two scalar functions X^iv) and H(k)(v) 



Sg^o 
Sg m n 



0, 



( X (k)(v)P m n + m^Qrn 1 ) + CC. 



(3) 
(4) 



where P m n and Q m " stand for scalar harmonics |2] , 

1 



p n 



-Q (fc) (x)<5 m ", 



1 



K 



O n 



(5) 
(6) 



<5(fe)(x) is the complex solution of the Helmholtz equation Q( fe )(x). TO ;m = — (A; 2 — 
K)Q(k)( x )> x = {X:S,f}, and the amplitudes X(k)( r l) an d A 1 (fe)( ? 7) satisfy the 
system of two second order equations convoluted to each other 



K 



/#,(•/) + [2 + 3*5(11)] ^{ k) (v) 



[\k)(v) + V(k)(v)] = o, 



+- 



k 2 - AK 



1 + 3c 2 M] [\ {k) ( v ) + m ( v )] = 0. 



(7) 



(8) 
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c s (ri) denotes the sound velocity: c s {rj) — \/ p 1 {rf) / e' {rf) 7^ 0. Then the density 
contrast 1 S = Se/e is a linear combination |5] 

s (v, x) = 5 (k) M Q(k) (x) + c.c, (9) 



and 



A 



(fe) 



3e(rj)a 2 (Tj) 



(/c 2 -4if)(A (fc) (r 7 )+ m (r 7 )) 



(10) 



where -4(fe) are arbitrary complex numbers. The system Q|SJ) defines the four 
dimensional phase-space, which means that the system possesses two physi- 
cal and two gauge degrees of freedom. Pure gauge solutions are known, and 
for A( fe ) (77) they respectively read [2] 



Gx{v) = 1, 

G 2 {rf) = ~{k 2 -K) 



a{rf) 



dij. 



(11) 
(12) 



Difficulties with fixing the gauge freedom have stimulated elaboration of the 
gauge-invariant theories |SJ E| |S1 |U] • After the famous Bardeen's paper JU| 
these theories were intensively developed for almost two decades. Yet, the orig- 
inal Lifshitz-Khalatnikov formalism, when appropriately extended, provides an 
equally good description of inhomogeneities. The procedure is as follows: we 
reduce (0-[8J) system to the single fourth-order equation for X/^ {rf) . 



+ { -5K + a 2 (r]) 



p{rf) )+{9H 2 {r ] )-e{r 1 ))c 2 s {v) 



+ (fc 2 - K) c 2 { v ) J X'( k) ( v ) + a{rf)H{rf) [-a 2 { v ) [e(r?) + 3p(v)] 

+ 2a 2 {r,)H 2 ( n ) [l + 3c 2 (77)] + (fc 2 - K) c 2 (r?)] A' (fe) (r?) = 0, (13) 

where H{rf) stands for the Hubble parameter H{rf) = a 1 {rf) / a 2 {rf) . The knowl- 
edge of the gauge solutions I|11H12|I enables one to extract the gauge space from 
the space of all solutions. First, we write the solutions in the form 

\k) {rf) = h Ai(fc) {rf) + / 2 A 2 ( fc ) {rf) + QiG 1(k) {rf) + g 2 G 2 (k) {rf) (14) 

with explicitly separated linear subspace giGi^{rf) + <72G 2 (fc)(?7) carrying all 
the gauge freedom (51,32 two arbitrary coefficients). Subsequently, we adopt 
the Darboux transform 2 I12| of {rf) to express the two remaining (physical) 

1 In more rigorous notation one should write 5 = Se/eo, where eo denotes background, 
unperturbed energy density. For the sake of simplicity we skip the subscript "0" in formulas. 

2 Darboux transform of a function f(x) is defined as f(x) = A(x)f(x) + B(x)f (x) (compare 
also formula 1181 in section 3) with arbitrary ^-dependent coefficients A and B. The same 
transform is extensively exploited in the soliton theory 
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degrees of freedom. A new perturbation variable B/f.\(rj) appears 



d \k)iv) V d G 2(fc)W 



a(r?) 



d ^G 2W (ry)' 



dr? V d7 / G i(/c)( ? 7)A dr ? G i(fe)(?/) 
On strength of (|13|l the variable B^(jj) satisfies the second order equation 
B'( k) {n) -[l- 3c 2 M] a(r))H( V ) B{ k) ( v ) 



(15) 



^+c 2 (r ] )\a 2 (rj)e(r } )-(k 2 -K)c 2 s (rj) 



(16) 



Equation Ultifl docs not contain gauge the coefficients, therefore, the variable 
B(k) (v) is gauge-invariant. Now, the function A(^) (77) can be found by the inverse 
Darboux transform of the two linearly independent solutions to equation Ijl6(l 3 



\ {k) (r 1 ) = h(K-k 2 ) 
1 



a(rj) 



B i(k) (v) 
a{rf) 



drj 



+MK-k J 



a(rj) 



B 2(k)(v) 



d?/ 



dv + 9iGi(k)(r]) +g 2 G'2(fe)(77). (17) 



On can insert \k)(ii) to Q) and solve it algebraically to obtain the second 
unknown function ji^{rf). As a result, by fixing four numbers {/1, /2, <7i, 32} one 
uniquely determines the metric correction 5g^ v . The metric correction is not 
gauge-invariant and cannot be such in any perturbation formalism. Quantities 
containing time integrals of B (k) (77) would exhibit similar properties, while those 
containing B^(i]) and derivatives are obviously gauge-invariant. 



3 Families of gauge- invariant variables in the 
synchronous reference system 

The metric perturbation can be written in the form (|17l) if the solutions of the 
equation (|16fl arc explicitly known. However, even in the case when these so- 
lutions are not known one may employ the equation (|16|) to "convert" some 
non-gauge-invariant perturbation variables into gauge-invariant ones. This par- 
ticularly refers to perturbation measures based on hydrodynamic quantities, 
since they involve the metric derivatives and are free of the metric elements 
themselves. The procedure is following: 

1) For the non-gauge-invariant perturbation variable X we take the linear com- 
bination of X and its time derivative X' (the Darboux transform) 



X(t?,x) = d(r/) 



d_ 



X(r),xL) + c 2 (r))X(r),x.) 



(18) 



3 The integration constants in the inverse Darboux transform can be set arbitrary, because 
the gauge freedom is already guaranteed by free choice of the coefficients gi and 32 ■ 
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which means that 



X(v,x) =^X( fc) (?7)Q( fc) (x) +c.c, 



(19) 



and 



x (k){v) = c\{rj) 



(20) 



2) We apply the equation to formally express the coefficients Xf k -s(r]) in the 
form 

X {k) ( V )=F 1 ( V )^B (k) ( V )+F 2 (r 1 )B {k) (i 1 )+F 3 (i 1 ) J (21) 

3) We look for c\{rj) and C2(ry) such that Xt k \{rj) contains solely function Bi k \{rj) 
and its derivatives but not integrals. (Integral over -q in f2"T|) introduces arbi- 
trary constant c/ k \, which depends on the wave vector k. The variable X(rj,x) 
(formula I|1H||) would contain then terms being arbitrary functions of the space 
coordinate x. In this way integration constant en,.-) restores one degree of the 
gauge freedom.) The relevant condition F^rj) — leads to a linear equation 
for 02(77), and leaves ci(?y) free. Therefore, the new gauge-invariant variable 
jf(?7,x) is actually a ci (^-dependent family of variables, where c\{rj) is an ar- 
bitrary function of time. Below we show some examples. 

3.1 The density contrast 

On strength of © , (ZHH) an d ifTTjl the density contrast 8^) (v) can be expressed 
by use of the gauge invariant function B^(rj). 



H{rj) 



K 



3 H( V ) [ B {k) {n) A 

2 e(j7) J a(i?) 



(22) 



The contrast 8t k \(r)) is not gauge- invariant, as it is well known |l(Jj . The gauge 
freedom is caused by the last term in (|22|) containing the integral of (77). 
Let us introduce a new variable — the "modified density contrast" 



With the aid of i|16|) one obtains 



d r\ 



s (k)(v) + c 2 (j7)<5 (fc) (77) 



5 {k) (r ) )=F 1 ( V )-^B {k) (r ] )+F2(r ) )B {k) (r 1 )+F 3 (r 1 ) J ^y^, 



(23) 



(24) 
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where the function ^3(77) reads 4 



By setting 



c 2 (t?) = -3— In 



d h / 1 e'fo) 
dr) \a(r]) e(rj) 

1 e'(rj) 



C2(r?) 



d?7 V«(»7) / ' 



(25) 



(26) 



one eliminates the integral oiB^ (77) (F 3 (rj) = 0), hence 5^ (rj) becomes a gauge 
independent quantity 



a(r?)e(?7) d?y 



a(r?)e(ry) 



e'(7j) 



(27) 



This variable is still defined up to the factor 01(77) being an arbitrary function 
of time. 

3.2 The expansion rate contrast 

Identically we proceed with inhomogeneities in the expansion rate 5 W(^,x). 

50(ji, x) = £ W (fc) (77) Q (fe) (x) + c.c, (28) 



where 



a 3 (77) 



a 2 (?7)(p(77)+e(77)) 2 



Hjrj) 
a 2 (77) 
3 
2 



A' 



a 2 (7/)(p(r/) +£(77)) 2 



2 (?(»?) + <»»)) a2(r?) 



B (k)(v) 



drj . 



We introduce £t?(m (77) defined as a linear combination 



= £1(77) 



d 77 



and eliminate the integral of B^)(r]). We finally obtain 



d7y v ' d?7 



a (7/) 



(29) 



(30) 



(31) 



4 -F\(»)) and F2(rj) are irrelevant to the problem discussed here, therefore we skip writing 
them explicitly. 

5 i? = vf'-a is the expansion 1131 . 
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3.3 Other variables 



By employing the same procedure one may define gauge-invariant variables 
which combine inhomogeneities in Hubble flow with inhomogeneities of the en- 
ergy density 

1 e'(r?) 



?(?7,x) = 01(77) 



M»7> x ) - o 



3H'( V ) 



5d{q, x) 



(32) 



It is clear that functions F(6, 5$,?, . . . , e, H, , . .) of the gauge-invariant variables 
{6, 5$,?, . . .} and the parameters of the background universe {e, H, . . .} are also 
gauge- invariant. None of these variables is preferred by the perturbation theory 
itself. The choice of the variable depends on its physical meaning and usefulness, 
its relation to other physical concepts and theories, and eventually on the ability 
to construct relevant observables. 



4 Perturbation equations in the canonical form 



In all the examples above we constructed the gauge- invariant variables which re- 
fer to the synchronous system of reference. The function C2(r)) was uniquely de- 
termined for each S(rj), *&(rj) and X(r)) by the demand of their gauge-invariance. 
The freedom to make an arbitrary choice of the time-dependent factor c± (rj) still 
remains, and therefore, we are now able to construct gauge-invariant variables 
with particular dynamical properties. Consider the variable 8(rj) with Cxirj) 
chosen as 



c l {ri) = a\rf)H i {n)- i 



We obtain the gauge-invariant variable 6 

? (fc) (r ? ) = a(r ? )F 2 (r ? )A- 
v ' drj 



(33) 



(34) 



which satisfies the propagation equation of the form 



a (17) c s (?7) 
with the function a (77) defined as 

a(»?) = a (v)\ 



6{ k) (v) + c 2 M(k 2 -K)5 (k) ( v ) = 0, (35) 



/ 1 e{jj)+p{rj) 
c s (7 ? ) m 2 (rj) 



(36) 



Now we introduce the new time parameter d£ = c s (rj)drj (acoustic conformal 
time |14|). After this reparametrization the propagation equation reduces to 
d'Alembert equation 

VV,%x) = (37) 

6 Please note the misprint in paper 1141 formula (7), already corrected in gr-qc/0302091. 
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in the Robertson- Walker space-time with coordinates {£, x, 8, cp} and the metric 
form 



g» v = a 2 (£) diag 



sin^^x) S in^(Vg X ) . 2 



(38) 



Therefore, the modified density contrast <5(£,x) propagates in the Robertson- 
Walker space-time with the scale factor a(£) in the same manner as minimally 
coupled scalar field in the Robertson- Walker space-time with the scale factor 
a(rj) |15|. All the classical results obtained in the field theory on the curved 
space-times apply to the density perturbations in the expanding universe. One 
can introduce the Lagrangian 

£ = \sTSj,v (39) 

for field <5(£,x) and derive the equation of motion l|37|) by use of the Euler- 
Lagrange equations 

Vu^^£-^C = 0. (40) 
d6 t „ 85 

An appropriate change of the perturbation variables |14j reduces the propa- 
gation equation of any gauge-invariant theory to the form of the equation (|16|) . 
By use of the procedure discussed above one can construct canonical variables, 
and as a consequence, each of these formalism can be expressed in the Lagrange- 
Hamilton language. 

The authors thank L.M. Sokolowski for reading the manuscript and for con- 
structive critics. 
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